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Abstract. We present a combined study of the thermal expansion and the thermal 
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known spin-state transition in LaCo03 is strongly affected by the exchange of the R ions 
due to their different ionic radii, i. e. chemical pressure. This can be monitored in detail by 
measurements of the thermal expansion, which is a highly sensitive probe for detecting spin- 
state transitions. The Co ions in the higher spin state act as additional scattering centers for 
phonons, therefore suppressing the phonon thermal conductivity. Based on the analysis of 
the interplay between spin-state transition and heat transport, we present a quantitative model 
of the thermal conductivity for the entire series. In PrCoC>3, an additional scattering effect is 
active at low temperatures. This effect arises from the crystal field splitting of the 4f multiplet, 
which allows for resonant scattering of phonons between the various 4/ levels. 
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1. Introduction 

Cobalt compounds are of particular interest due to the possibility that Co ions can exhibit 
different spin states and hence the occurrence of temperature-driven spin- state transitions. The 
most prominent example is LaCo03, which has been intensively studied and controversially 
debated for more than fifty years, see e. g. EE 13 BIS El 171 El El- The Co 3+ ions in LaCo0 3 
feature a 3d 6 configuration which in principle can occur in three different spin states: a 
nonmagnetic low-spin (LS) (t% g e° g , S = 0), an intermediate- spin (IS) (1% e g , S = I) and a 
high-spin (HS) state (4 g e g' $ = 2). It is generally agreed that the Co 3+ ions in LaCo0 3 
realize the LS state at low temperatures. Above approximately 25 K a higher spin state, 
either IS or HS, becomes thermally populated affecting various physical properties, e. g. the 
magnetic susceptibility^ or the thermal expansion a, which both exhibit pronounced maxima 
in their temperature dependencies [fTOl ITTTl . The susceptibility is obviously affected because 
the excited spin state, either IS or HS, induces a strong increase of the magnetization above 
25 K. The thermal expansion is affected due to the different ionic radii of the smaller LS-Co 3+ 
with empty and the larger IS- or HS-Co 3+ ions with partially filled e g orbitals. The spin-state 
transition can be well described in a LS - IS scenario, i. e. the excited spin state is the IS state, 
with a constant energy gap of A Co = 185 K [[101 [TTJ. However, more recent investigations 
show, that a LS/HS model including spin-orbit coupling is more reasonable |[T2l [T3l [T4| . A 
consequence of the latter model is a temperature-dependent energy gap A Co (T), which strongly 
increases with increasing temperature lfT3ll . 

The spin-state transition in LaCo03 is strongly affected by both, heterovalent and 
isovalent doping on the La site. The former possibility causes hole doping and chemical 
pressure and suppresses the spin-state transition due to the implementation of Co 4+ ions 
[fT51 [T6l . The latter one, i. e. chemical pressure without changing the Co valence, is usually 
realized by introducing trivalent rare-earth ions R 3+ . In 7?Co0 3 the spin-state transition is not 
suppressed but its onset is shifted to higher temperature. The energy gap between the LS and 
the excited spin state increases from about A Co = 185 K for R = La to > 2000 K for R = Eu 
ifTTfl. Moreover, due to the decreasing ionic radius of the lanthanide series, the structure 
changes from rhombohedral in LaCo0 3 to orthorhombic for R = Pr, Nd and Eu. Recently, 
it has been reported, that the low-temperature thermal conductivity of LaCo03 is also very 
anomalous lfT71[T8Tl . This behaviour has been qualitatively attributed to the onset of the spin- 
state transition. However, a quantitative analysis of the anomalous thermal conductivity of 
LaCo0 3 has not been presented yet. 

The aim of this paper is to develop a consistent picture of the influence of the spin-state 
transition on the thermal conductivity. Therefore, we measured the thermal conductivity k(T) 
of the series 7?Co0 3 with R = La, Pr, Nd and Eu. Moreover, we studied the thermal expansion 
a(T), which is a very sensitive probe to investigate spin-state and also crystal-field transitions 
and their coupling to the lattice. In the quantitative analysis of our data, we will consider 
both models of the spin-state transition, which are favoured for LaCo0 3 in the literature, i. e. 
we will consider the LS-IS scenario with a constant A Co and the more recently proposed 
spin-orbit coupled HS (SOcHS) model with a temperature dependent energy gap A Co (T). 
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Table 1. Characteristic properties of the investigated crystals. Here, T on denotes the onset 
temperature of the spin-state transition, Ac the energy gap between the LS ground state and 
the excited spin state (either IS or HS) and y is related to their ionic radii difference (see text). 
Due to the large T on of EUC0O3 we can only give lower limits for the values of y and Ac , as 
it was also the case in the related analysis of the magnetic susceptibility (TTJ. The last three 
columns give the room-temperature values of the thermopower Srt, which for the LaCo0 3 
crystals have been taken from lfT8l . and the crystal dimensions (sample cross section A and 
sample length L) which are important for the measurements of k (the cross section of S5 is an 
approximate value, since it is not of rectangular shape). 



Sample 


T on (K) 


Aco(K) 


y(%) 


S R tOuV/K) 


A (mm 2 ) 


L(mm) 


LaCo0 3 (SI) 


25 


185 


0.7 


-700 


0.8 x 1.5 


3.7 


LaCo0 3 (S2) 








1000 


1.1 x0.8 


2.3 


LaCo0 3 (S3) 








-600 


1.4x0.8 


3.4 


LaCoC-3 (S4) 








1000 


0.8 x 1.4 


2.6 


LaCo0 3 (S5) 








-300 


=* 2.3 


4.7 


PrCo0 3 


175 


1200 


2.8 


-400 


1.0x3.8 


2.5 


NdCo0 3 


230 


1700 


4.8 


-400 


1.7 x 2.0 


2.8 


EuCoQ 3 


400 


> 2000 


> 10 


-500 


0.3x0.85 


2.0 



2. Experiment 

All 7?Co0 3 crystals have been grown in a floating-zone image furnace. We examine five 
different LaCo0 3 crystals identical to those used in [18], where details of the sample 
preparation and characterization are given. For EUC0O3 this information can be found in 
IfTTI . The NdCo0 3 and PrCo0 3 single crystals have been grown in the same way as those 
of LaCo0 3 and EuCo0 3 . Characteristic properties of all crystals are listed in table [TJ The 
thermal conductivity measurements have been performed by a standard steady-state method 
using a differential Chromel-Au+0.07 % Fe-thermocouple [fT9ll . The thermal expansion below 
« 200 K was measured using a home-built high-resolution capacitance dilatometer [|20l . 
whereas the high-temperature measurements 135 K < T < 670 K were performed using 
commercial inductive dilatometer (TMA 7, Perkin-Elmer). 

3. Results and Discussion 

3.1. Thermal Expansion 

In figure[l](a) we show the thermal-expansion coefficients a(T) = l/L-dL/dT of i?Co0 3 . The 
low-temperature results of LaCo0 3 and EuCo0 3 were already discussed in detail in [TTOlfTTTl . 
The thermal expansion consists of a phononic part and a Schottky contribution, caused by 
the thermal population of the e g orbitals of the Co 3+ ions, i.e. a = arp h + ars C h- The latter 
contribution causes the large maximum at 50 K in LaCo0 3 . To further analyze the data 
we subtract orph, which we estimate using a Debye function with the Debye temperature 
©d = 600 K of LaCo0 3 j[2T|. Here, we used the same D for the various cobaltates, 
but sample-dependent prefactors determined by scaling the Debye function to the low- 
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T (K) T (K) 

Figure 1. (a) Thermal expansion of RC0O3 with R = La, Pr, Nd and Eu. The solid line is 
the estimated phonon contribution of NdCoC>3 (see text), (b) Anomalous part of the thermal 
expansion obtained by subtracting the phononic background. For clarity the different data 
sets are offset by 10~ 5 /K with respect to each other. The arrows signal the approximate 
onset temperature of the spin-state transition and the solid circles denote the metal-insulator 
transition temperature estimated from resistivity measurements for each compound. The solid 
lines are fits of the respective Schottky anomalies, see equations (fl]i and (O, assuming a 
constant energy gap Ac between the different spin states of Co 3+ . The dashed line for LaCoC>3 
is a similar fit using Ac (T) obtained from the magnetic susceptibility (see text). 



temperature data of each compound. As an example, we show a P h for NdCo03 in figure[H(a); 
arph of EUC0O3 (PrCoOs; not shown) is slightly larger (smaller). Since a clear separation of 
Qfph and ar Sc h is not possible for the low-temperature data of LaCo0 3 , we used ar Ph of EuCo0 3 
also for LaCo03. The resulting as C h of all crystals are shown in figure [H(b). We note that, in 
particular for PrCo0 3 and NdCo0 3 the spin-state transition is seen much better in the thermal- 
expansion data than in the magnetic susceptibility. The reason is that a^h is rather small 
compared to the total thermal expansion, whereas x is dominated by the large contribution of 
the Af moments of the Pr 3+ and Nd 3+ ions, respectively, which makes a further analysis rather 
uncertain. The insulator-metal transitions occurring above about 500 K, see e. g. |fTTll22)| also 
cause anomalies in a(T). The solid circles in figure [H(b) signal the transition temperatures 
T M i, which have been determined from resistivity measurements on our crystals (not shown). 
Within the LS/IS scenario with a constant A Co , a Sc h(T) is given by IflOl 

Aco vexp(-A Co /r) 
T 2 (l+vexp(-A Co /r)) 2 
where v = 3 is the degeneracy of the excited spin state and y is a measure of the ionic radii 
difference of the Co 3+ in the LS and the excited spin state. Moreover, y is related to the 
uniaxial pressure dependence of A Co via y = k B /Vf u ■ dA Co /dp a , where k B is Boltzmann's 
constant, V fu ^ 56 A 3 the volume per formula unit and p a means uniaxial pressure along the 
direction of which a is measured [|23ll24l . To get the hydrostatic pressure dependence one has 
to use the volume expansion which in the case of the twinned 7?Co0 3 single crystals is given 
by 3 a. The corresponding fit of a of LaCo0 3 (already presented in IfTOl ) is shown by the solid 



«sch(D = 7 — — — —7 ' (!) 
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line in figure [U(b). The fit parameters are A Co = 185 K and y = 0.007 giving a hydrostatic 
pressure dependence d(ln A Co )/<9/? h ydr - 45 %/ GPa, in agreement with the increase of A Co of 
42 % obtained from measurements of^ under ^1.1 GPa [|25l . 

For the SOcHS model one has to consider the temperature dependence of A Co ffT3ll . which 
can be calculated from the measured x(T) [11]. For LaCo03, this yields a linear increase 
AS^CT) = A9, + a • T with A° o = 135 K and a = 1.66 in the temperature range almost up to 
the MI transition. Using this AS^(r) in the partition sum, equation (QQ) modifies to 



The corresponding fit of asch with y s as the only free parameter gives the dashed line in 
figure Q](b). Obviously, both fits hardly differ because the modified A^ is compensated by a 
larger value of y HS = 0.02. We note, however, that this does not necessarily correspond to a 
larger pressure dependence of A9, , because ar Sc h of equation © is determined by the pressure 
dependencies of A° o and that of the slope a; see e. g. the discussions in [|23ll26ll27l . In view 
of the small differences between both fits for LaCo0 3 , and since there are no indications for a 
temperature-dependent A Co from other physical quantities, the fits of asch °f me other i?Co0 3 
crystals have been done for constant energy gaps only. Since equations (0Q) and © are identical 
for a = 0, this analysis is not able to distinguish between both models. The obtained values of 
A Co and y are given in table CD As expected, A Co strongly increases with decreasing radius of 
the R 3+ ions and therefore the spin-state transition monotonically shifts to higher temperature. 
The onset temperature r on of the spin-state transition for each compound, T on 25, 175, 230 
and 400 K for R = La, Pr, Nd and Eu, respectively, are marked by arrows in figure [TJ(b). 

3.2. Thermal conductivity 

Figure [2](a) displays the thermal-conductivity data of i?Co0 3 . As has been already found in 
previous studies on LaCo0 3 ifTTl [T8ll . the overall shape and qualitative behaviour of k{T) is 
rather unusual. First of all, k is rather low in the whole temperature range and its temperature 
dependence clearly deviates from the typical behaviour expected for a conventional phononic 
heat conductor. The thermal conductivity of LaCo0 3 exhibits a maximum around 20 - 30 K. 
Towards higher temperatures k rapidly drops and features a minimum around 150 K instead 
of the expected 1 /T-decrease (see also |[T8l0 . Qualitatively, this minimum can be traced back 
to the spin-state transition of LaCo0 3 , which sets in close to the maximum of k. The thermal 
population of the Co e g levels induces a certain fraction of Co 3+ ions with a larger ionic radius 
compared to the LS Co 3+ ions. These randomly distributed larger Co 3+ ions do not only cause 
the huge anomaly in a(T), but also lead to additional disorder in the lattice and therefore act 
as additional scattering centers for phonons. Hence, the thermal conductivity is additionally 
suppressed. A quantitative description will be given below. 

Since the spin-state transition shifts to higher temperature when La 3+ is replaced by 
smaller R 3+ ions, the above scenario suggests that the onset of the suppression of k(T) of 
i?Co0 3 systematically shifts towards higher temperature with decreasing radius of the R 3+ 



asch(T) = y 



HS 





vexp(-(A« o + a7/)/r) 



(2) 




Figure 2. (a) Thermal conductivity of RC0O3 with R — La, Pr, Nd and Eu; the solid line is a 
fit via equation (01. (b) Thermal conductivity of five different LaCoC>3 single crystals SI -S5 
(note the different scale). The inset shows the maximum values of the thermal conductivity 
Kmax of these crystals vs. the saturation values of the magnetization M sat measured at 2 K. As 
shown by the sold line, we find an approximate proportionality K max oc 1 /M sat (see text). 



ion. Let us first consider the data of EuCo0 3 , which is the compound with the smallest R 3+ 
ion and the largest Ac ^ 2000 K [11 J. Hence, the spin- state transition should not affect 
the thermal conductivity below room temperature and, indeed, the thermal conductivity of 
EuCo0 3 exhibits much larger values than in LaCo0 3 . Its low-temperature maximum exceeds 
that of LaCo0 3 by about one order of magnitude. Moreover, the temperature dependence of 
k of EuCo0 3 meets the expectation of a purely phononic heat conductor. 

The thermal conductivity of NdCo0 3 is lower than k of EuCo0 3 in the entire temperature 
range and, at a first glance, its overall temperature dependence seems to be consistent with a 
conventional phononic picture, too. A closer inspection of the data yields, however, a slight 
slope change towards a somewhat steeper decrease of k(T) at « 230 K, i. e. close to the onset 
temperature of the spin-state transition, see figure[H(b). This slope change supports the above 
proposed explanation that the spin-state transition causes additional scattering centers. 

PrCo0 3 also exhibits an unusual temperature dependence of k, which appears somewhat 
similar to that of LaCo0 3 . Above the low-temperature maximum, k(T) rapidly decreases to a 
small value of ^ 3.5 W/Km and remains almost temperature independent between about 80 K 
and 175 K. Above this temperature a further decrease of k sets in. Since this decrease starts 
again at the onset of the spin- state transition (see figure [H(b)), it is natural to attribute it to the 
additional scattering due to the presence of Co 3+ ions in different spin states. In contrast, the 
anomalous temperature dependence of k(T) at lower temperature is not related to the spin- 
state transition and, as will be discussed separately in section H31 this anomalous feature is 
related to the crystal-field splitting of the Af shell of Pr 3+ . 

Summarizing the presented data so far, we observe in all crystals with r on < 300 K that 
this onset coincides with an additional suppression of k. Since Ac strongly increases with 
decreasing R 3+ radius lfTTll2"2"l . T 0B also shifts towards higher temperature and the additional 
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suppression of the thermal conductivity is less obvious, because at higher temperature the 
usual suppression of k due to phonon-phonon Umklapp scattering becomes more and more 
dominant. However, the data in figure [21(a) also show, that k is not only suppressed at 
temperatures above T on . Instead, with increasing size of R 3+ a monotonic decrease of k in 
the entire temperature range is found (here we have neglected the additional suppression of k 
in PrCoC^). If the spin-state transition would be the only effect causing additional disorder 
scattering, k(T) should approach the 'bare' phononic value below T on . Taking EuCo0 3 with 
stable LS Co 3+ ions up to above room temperature as a rough measure for this 'bare' phononic 
heat conduction, it becomes clear that this expectation is not observed in our data. 

In general, the magnitude of the low-temperature maximum of k depends on the sample 
quality, which determines the maximum free path. Thus, the different values within the 
7?Co03 series could arise from different sample qualities. However, there is no indication 
for a systematic increase of the crystal quality with decreasing R 3+ ion radius @ In order to 
study a possible influence of the crystal quality, we measured k of five different LaCo0 3 single 
crystals. As shown in figure 12(b) the high-temperature behaviour is essentially the same for 
all crystals, but the low-temperature maxima differ considerably. Nevertheless, all of them 
remain much lower than that of EUC0O3. We relate the low-temperature maxima of k to the 
defect concentrations by considering the low-temperature magnetization, since a really pure 
LaCo0 3 with all Co 3+ ions in the LS state would be non-magnetic (apart from core dia- and 
van Vleck paramagnetism). At 2K we measured magnetization up to 14 Tesla in order 
to reach its saturation M sat , which is a measure of the content of magnetic impurities. The 
inset of figure [2] displays the low-temperature maxima of k(T) vs. M sat , which indeed roughly 
follow a 1 /M sat dependence as shown by the solid line. 

The most likely source for magnetic impurities in LaCoC>3 is a weak oxygen 
nonstoichiometry, which changes the average valence from Co 3+ towards Co 2+ or Co 4+ . 
Due to the odd number of 3d electrons, the latter are magnetic in all possible electronic 
configurations, but in most cases Co 2+ and Co 4+ realize t 5 lg e 2 g and t 5 2g e° g states with spin 5=3/2 
and S = 1/2, respectively. Experimental evidence for a weak charge-carrier doping in 
nominally pure LaCo03 stems from the large values of different signs of the room-temperature 
thermopower for the different LaCo0 3 crystals; cf. table [Hand the discussion in lPT8ll . From 
magnetization data on Sr-doped LaCo03 it is known that low doping concentrations induce 
so-called high-spin polarons with effective spin values up to S = 16 Il28l . Such polarons 
are formed around the doping-induced Co 4+ ions, which cause a spin-state transition in the 
adjacent Co 3+ ions. Most probably, this is a transition from the LS to an IS state, because 
the electronic configuration t 2g e l g of Co 3+ is favourable for the hopping of the e g electron for 
both kinds of neighbours, Co 4+ with t\ g e° g and Co 2+ with t\ g e 2 g . Other possible configurations 
often only allow a t 2g hopping or may cause the so-called spin-blockade effect (291 EDI. The 

§ One may speculate that the low-temperature maximum of k is also influenced by twinning. All these single 
crystals are twinned, but it seems plausible that the average size of the twin domains increases from LaCo03 to 
EuCoOi, because the deviation from the cubic symmetry increases. However, the observed correlation between 
Kmax and the magnetic impurity concentration for the different LaCo03 crystals, makes it rather unlikely to 
explain the different low-temperature maxima of k in terms of different domain sizes. 
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formation of high- spin polarons will depend on A Co between the LS ground state and the 
excited spin state of Co 3+ . Since A Co strongly increases with decreasing R 3+ radius, the 
occurrence of such polarons in the i?Co0 3 series should get more and more unlikely for 
smaller R 3+ ions. Unfortunately, the presence of high-spin polarons in PrCo0 3 and NdCo0 3 
cannot be studied by magnetization measurements because of the strong magnetism of the 
partially filled Af shells of Pr 3+ and Nd 3+ . For EuCo0 3 such an analysis is again easier, 
because Eu 3+ only exhibits van Vleck paramagnetism and, in fact, the impurity contribution to 
the magnetic susceptibility in EuCo0 3 is one order of magnitude smaller than in LaCo0 3 |jTT||. 
Moreover, the magnetic impurity contribution of the Lai_ A Eu T Co0 3 series systematically 
decreases with increasing x, i. e. with decreasing average ionic size on the La site and with 
increasing A Co [fTTTl . These observations strongly support the above idea that the systematic 
dependence of the low-temperature thermal conductivity of i?Co0 3 on the ionic size of R 3+ 
is related to phonon scattering by high-spin polarons, whose formation depend on a weak 
oxygen-nonstoichiometry and on A Co . 

4. Quantitative analysis of the thermal conductivity 

For a quantitative analysis of the thermal conductivity we use an extended Debye model OTI 

A 1 &D/T 

k T C x e x 



Here, D denotes the Debye temperature, v s the sound velocity, co the phonon frequency, 
x = hco/k^T and r(x, T) the phonon relaxation time. The scattering rates of different 
independent scattering mechanisms sum up to a total scattering rate 

t~\x, T) = T~ h l + Tpt 1 + r;i = V -^+Pco 4 + UTco 3 exp(%). (4) 

The three terms on the right-hand side refer to the typical scattering rates for phonon heat 
transport in insulators, namely boundary (r b d; L = 1 mm is the characteristic length scale of 
the sample [32] and was kept fixed for all crystals), point defect (r pt ) and phonon-phonon 
Umklapp scattering (r um ). Since the lattice spacing is a natural lower limit for the mean-free 
path, k{T) is expected to approach a minimum value at higher temperatures. In order to model 
this, t(x, T) in equation © is replaced by max{r-z(x, T), £ m i n /v s } with the minimum mean-free 
path £ m [ n being of the order of the lattice constant; compare the discussion in 11331 . 

As shown by the solid line in figure El the thermal conductivity data of EuCo0 3 can 
be reproduced reasonably well by equation©. In order to reduce the number of adjustable 
parameters, we used D = 600 K and v s = 3900 m/s determined for LaCo0 3 [|2T1 [34l . The 
other parameters have been adjusted to fit the experimental data. We obtain P = 3.1 ■ 10~ 43 s 3 , 
U = 4.9 • 10" 31 s 2 /K and u = 7.6; € mia is not reached in EuCo0 3 up to 300 K. 

4.1. Scattering related to the spin-state transition 

In the following analysis we consider the thermal conductivity of EuCo0 3 as the bare 
phononic heat conductivity that is not influenced by the spin-state transition. For the other 
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members of the 7?Co0 3 series we introduce an additional scattering mechanism arising from 
the thermal occupation of the Co 3+ e g shells, which causes additional lattice disorder due to 
the presence of smaller and larger Co 3+ ions in the LS and the excited spin state, respectively. 
Considering the Co 3+ ions of different sizes as point-like defects, we assume an oj a frequency 
dependence for this additional scattering rate (cf. equation ©) 

r^ = C-f(T)-co\ (5) 

where C is an adjustable parameter describing the scattering strength and f(T) models the 
temperature dependence of the additional lattice disorder. 

The unknown functions C • f(T) for LaCo0 3 , PrCo0 3 and NdCo0 3 can be calculated 
from the experimental data in the following way: we use equation © with all the parameters 
obtained from the fit of k of EuCo0 3 kept fix and add r"j^. Then the total r is used in 
equation © and the values of C • f(T) at each temperature are determined point by point 
from the measured k(T). The resulting C • f(T) curves for all three i?Co0 3 samples are shown 
by the open symbols in figure [H In order to avoid additional complications arising from the 
aforementioned high-spin polarons, we restrict the analysis to the temperature range above 
25 K. For LaCo0 3 the resulting C • f(T) has a pronounced maximum around 100 K, and 
strongly resembles the temperature dependence of the magnetic susceptibility of LaCo0 3 . 
Here, we have used Sample S2 which exhibits the highest value of the low-temperature 
maximum of k, i. e. the sample which is least affected by impurity scattering. We note, 
however, that in the considered temperature range above 25 K the resulting C • f(T) is not 
very different for the other LaCo0 3 crystals. For PrCo0 3 , C • f(T) features a clear maximum 
followed by a slope change around 175 K. The obtained values of C • f(T) for NdCo0 3 are 
much lower than those of LaCo0 3 and PrCo0 3 and there is no maximum, but again a clear 
slope change around T on 230 K. Thus, we find in all three RCo0 3 compounds either an 
increase of C • f(T) or a slope change towards larger values at temperatures close to the 
respective onset temperatures of the spin-state transition; see figured] 

In order to model the additional scattering rate we use the so-called Nordheim rule 113511 . 
which is a successful approach to describe k in disordered mixed alloys. There, the additional 
scattering rate is proportional to x(l - x), with x being the fraction of one of the constituents. 
In the present case this corresponds to the fractions rci S / H s of Co 3+ ions in the higher and 
n Ls = (1 _ "is/hs) m the LS state. Thus, the temperature dependence of the additional scattering 
rate of equation © is expected to be given by (see equations ® and ©) 
= vexp(-A Co( r)/r) 
(l + vexp(-A Co (r)/r)) 2 

The solid line in figure [3](a) corresponds to a calculation of C • f(T) via equation © 
using the temperature-independent Ac = 185 K of the LS-IS scenario. Here, C is the 
only adjustable parameter, and for C = 5.9 • 10 _40 s 3 this model yields a good description 
of the experimental data in the temperature range from 30 K to 150 K. The dashed line shows 
the corresponding fit (yielding C = 11 • 10" 40 s 3 ) for the same temperature range based on 
the SOcHS scenario with a temperature-dependent A Co (T), which is in reasonable agreement 
with the data, too. Obviously, both models are not able to describe the decrease of f(T) above 
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T (K) T (K) 

Figure 3. Open symbols denote the experimentally obtained additional scattering rates 
C ■ f(T), which are needed to suppress k from a hypothetical bare phonon heat conductivity 
represented by k of EUC0O3 to the measured values of k of LaCoC>3, PrCo03 and NdCoC>3(see 
text for details). In all three RCoOt, an additional scattering sets in above the respective 
onset temperature T oa (marked by the arrow in each panel) of the spin-state transition. The 
lines show the expected additional scattering rates due to the spin-state transition for a LS/IS 
scenario (solid) and for the SOcHS model (dashed); see equations (O and ([6). In panel (a) the 
scattering strengths C for each model have been adjusted by fitting the lines to the experimental 
data in the temperature range below 150 K. These values have been kept fixed to calculate the 
lines in panels (b) and (c), and, as shown by the solid symbols, the additional scattering above 
T on completely vanishes, when these lines are subtracted from the experimental C ■ f(T). 

150 K. This disagreement does, however, not contradict the models, because for T > 150 K the 
measured k of LaCo0 3 is already close to its minimum value arising from the aforementioned 
lower limit of the mean-free path. For C • f(T) of PrCo03 and NdCo03, we keep the values 
of C fixed and use the gap values A Co = 1200 K and A Co = 1700 K obtained from the analysis 
of a(T), cf. table [Q The corresponding curves, which for the different models only differ 
in the scattering strengths, are shown as solid and dashed lines in Figs.[3](b) and (c). When 
we subtract these lines from the experimental data, the slope changes around r on completely 
vanish, as is shown by the solid symbols in Figs.[3](b) and (c). This gives further evidence that 
the spin-state transition is indeed responsible for the suppression of the thermal conductivity 
of PrCo03 and NdCo03 as it is the case for LaCo03. 

Encouraged by this finding we tried to fit the thermal conductivity of LaCo0 3 via 
equation© including the additional scattering rate given by equations © and © in 
equation ©. Again, we restrict our analysis to T > 25 K and keep most of the fit parameters 
fixed to their values as obtained from the Debye fit of the bare phononic k of EuCo0 3 . 
From these parameters, we only readjust the prefactor P of the point-defect scattering term 
in equation ©, because the number of point defects varies in different crystals, and point- 
defect scattering also reduces k at high temperatures. In addition, the scattering strength C is 
readjusted, but changes only little (to C = 6 • 10~ 40 s 3 and 12 • 10~ 40 s 3 for the LS/IS and the 
SOcHS model, respectively), and we consider the minimum mean-free path £ min . The results 
for the LS - IS as well as the SOcHS model are displayed in figure IU Both models describe 
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Figure 4. Thermal conductivity (symbols) of LaCoC>3 with fits based on equation (0 assuming 
(a) the LS - IS or (b) the SOcHS model for the spin-state transition. The dashed (solid) lines 
are obtained when a minimum mean-free ( m i n is (is not) taken into account (see text). 



the thermal conductivity well over a large temperature range of about 200 K, in particular, 
when the minimum mean-free path is taken into account. As was already seen in figure |3] 
the quality of the fits based on the different models hardly differs. Therefore, the analysis of 
the thermal conductivity of i?Co0 3 does not allow to judge, whether the LS/IS or the SOcHS 
scenario is preferable. 

4.2. Scattering related to crystal-field excitations 

Next we consider the additional suppression of k of PrCo03 below r on , which is not related to 
the spin-state transition. As can be seen in figure [H the thermal expansion of this compound 
also exhibits an additional anomaly at lower temperatures. This anomaly arises from the 
crystal-field splitting of the Af shell of the Pr 3+ ions. In an orthorhombic crystal field the 
3 H 4 multiplet of the free Pr 3+ ion splits into nine singlets. We are not aware of experimental 
investigations of the crystal-field splitting in PrCo03, but it is known that the crystal-field 
levels of other orthorhombic PrA0 3 are very similar Il3~6l [371 [38ll . Only the lower levels 
slightly depend on the exact orthorhombic distortion, which is mainly characterized by the 
A-O-A canting angle. The energy splitting between the crystal-field ground state and the 
first excited level in PrCo0 3 can be determined from a fit of the low-temperature thermal 
expansion by a Schottky anomaly. As shown in figure 0(a), such a fit for two non-degenerate 
levels (setting v = 1 and replacing A Co by A CF in equation (OQ)) describes the experimental data 
well and yields Acf = E\ = 70 K and y = 0.056 %. For the higher- lying energy levels we 
use the values E 2 = 151 K, £ 3 = 174 K, E 3 = 235 K, . . . of PrNi0 3 (381, which has almost the 
same A-O-A canting angle as PrCo0 3 . We note here that the thermal population of these 
higher-lying levels may contribute to a(T) at higher temperatures, but below about 50 K the 
two-level system used is sufficient to model the experimental data. 

The most likely source for the suppression of the thermal conductivity of PrCo0 3 in 
the intermediate temperature range is resonant scattering between different 4f crystal-field 
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Figure 5. Low-temperature (a) thermal expansion and (b) thermal conductivity of PrCoC>3. 
The thermal-expansion data (symbols) in panel (a) are well described by a Schottky anomaly 
(line) of a two-level system with non-degenerate levels with energy splitting E\ - Eo = 70 K. 
The lines F01, F02 and F12 in panel (b) are fits for T > 20 K to the thermal conductivity data 
(symbols) including either the resonant scattering process Eo «-» E\, Eo <-» £2 or E\ <-> E2, 
respectively. A sketch of the energy levels is also included 



levels, as it is also known from e.g. rare-earth garnets ll39ll . In such a resonant scattering 
process a phonon is absorbed by inducing a transition between different energy levels and 
then re-emitted in an arbitrary direction (see e. g. 11401 HTTP , causing an additional thermal 
resistance. A mechanism based on random volume changes as discussed above for the spin- 
state transition of Co 3+ is unlikely because the Af orbitals are inner shells. This difference is 
directly reflected in the magnitudes of the Schottky peaks, which are mainly determined by 
the underlying volume change. The low-temperature peak in a of PrCo0 3 is e. g. one order of 
magnitude smaller than that of LaCo0 3 ; see figured] and H] 

In order to model k(T) in the low-T range we include the rate 

r- = R- 2-j • (Ni + Nj) (7) 

(A^ F -- 2 ) 

in equation© describing a direct resonant scattering process Il40ll42l . Here, A C f = \Ei - Ej\ 
is the energy difference between the two levels, and Ni, Nj denote their population determined 
by the partition sum. For the fit we keep again D = 600 K and v s = 3900 m/s, while P, U, u 
and R are adjusted with respect to the data. The fits have been restricted to T > 20 K, since 
other processes may become dominant at very low temperatures. 

In a first attempt, resonant scattering between E and the first excited level Ex is 
considered, yielding the fit labelled F01 in figure |5](b). This fit only reproduces the data 
for T > 100 K but totally misses the low-temperature maximum in k(T). Basically the same 
result is achieved by considering the process E «-» E 2 (curve F02 in figure 13(b)). A much 
better description of the data is obtained, however, for the process E\ «-» E 2 (curve F12; with 
parameters P = 5 • 10~ 43 s 3 , U = 4.85 ■ 10- 18 s 2 /K, u = 5 and R = 3.8 • 10" 40 s 3 ). Hence, the 

|| We also mention that our attempts to describe the additional damping of k due to the spin-state transition of 
LaCoC>3 by a resonant scattering process completely failed. 
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transition <-> E 2 appears to be the dominant scattering channel of the phonons in PrCo0 3 
at intermediate temperatures. 

5. Summary 

In this paper we have presented measurements of the thermal-expansion coefficient and the 
thermal conductivity of 7?Co0 3 , with trivalent R = La, Pr, Nd and Eu. As a consequence 
of the chemical pressure due to the decreasing ionic radius from R = La to Eu, the energy 
gap between the LS ground state and the first excited spin state of the Co 3+ ions increases 
drastically and the temperature-induced spin-state transition of the Co 3+ ions is systematically 
shifted towards higher temperature. This strong pressure dependence is also reflected in 
pronounced Schottky anomalies in the thermal-expansion coefficients, making this quantity 
a sensitive probe to analyze spin-state transitions. The microscopic origin of the pressure 
dependence is the volume change due to the partial occupation of the e g levels of the Co 3+ 
ions in the excited spin state in contrast to the empty e g levels of LS Co 3+ ions. Another 
consequence of this volume difference is a strong suppression of the phonon heat transport, 
which sets in above the onset temperature of the spin-state transition, because the random 
occupation of the Co sites with Co 3+ ions in different spin states, i. e. of different sizes, causes 
additional lattice disorder. By considering this disorder as point-like defects, our analysis 
clearly shows this effect in the temperature dependence of the corresponding scattering 
rates of LaCo03, PrCo03 and NdCo03, where the spin- state transition sets in below room 
temperature, in contrast to EuCo0 3 which exhibits a stable LS state up to about 400 K. 

Quantitatively, the temperature dependence of the additional scattering due to the spin- 
state transition can be described in a model based on the Nordheim rule, which is usually 
applied to describe scattering in disordered mixed alloys. For the cobaltates, this means that 
the temperature dependence of this additional scattering rate is given by the product of the 
occupation numbers of the LS and of the excited spin states. Moreover, we find that the 
scattering strength is constant for LaCo0 3 , PrCo0 3 and NdCo0 3 . Our analysis does, however, 
not allow to distinguish between the different models discussed for the spin-state transition in 
LaCo0 3 , namely a LS-IS scenario with a temperature-independent energy gap or the more 
recently proposed model based on a spin-orbit coupled high-spin state with a temperature- 
dependent energy gap. 

The above described model explains a damping of the phonon heat transport above 
the onset temperature of the spin-state transition. Experimentally, however, we find that 
as a function of increasing size of the R 3+ ions the thermal conductivity decreases also in 
the temperature range below the respective onset of the spin-state transition. By studying 
five different LaCo0 3 single crystals, we could show that the low-temperature thermal 
conductivity inversely scales with the amount of magnetic impurities. Most likely, these 
impurities arise from a weak oxygen nonstoichiometry giving rise to a certain amount of 
Co 2+ or Co 4+ ions. Since it is known that these magnetic Co ions tend to form so-called 
high-spin polarons by inducing a spin-state transition in the neighbouring Co 3+ ions, our data 
suggest that the phonon damping at low temperature arises from such a doping-induced spin- 
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state transition. The drastic increase of the energy gap to the excited spin state of Co 3+ for 
decreasing R 3+ ionic radii makes the formation of spin polarons more and more unlikely, 
and thus can naturally explain the systematic increase of the thermal conductivity in the 
temperature range below the temperature-induced spin-state transition. 

Additional low-temperature features have to be considered in PrCo03, because the 
thermal expansion exhibits another Schottky anomaly and there is an additional suppression 
of the thermal conductivity. Both effects are not related to the temperature-dependent spin- 
state transition of Co 3+ , but stem from the crystal-field splitting of the ground-state multiplet 
of the Af shell of Pr 3+ . Our analysis demonstrates that the Schottky anomaly of the thermal 
expansion arises from a thermal occupation of the first excited crystal-field level, whereas the 
phonon damping is a consequence of resonant phonon scattering between the first two excited 
crystal-field levels. 
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